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-1- SDS/1965 
Radar Tracking System 
This invention relates to radar tracking 
systems for missiles directed against airborne targets. 

The radar system may operate in an active mode 
in which the target is illuminated by radiation from the 
5 missile and in which reflected radiation from the target 
(termed "skin echo") is received by the missile during 
its flight. The reflected radiation is processed to 
obtain the bearings , azimuth and elevation, and the 
speed or range of the target so that the missile can 
10 follow the changes of direction and speed of the target. 

An object of the present invention is the 
improvement of terminal accuracy in a pulse doppler 
radar seeker. 

According to one aspect of the present 
15 invention, in a monopulse radar tracking system 

including either or both of a doppler tracking loop and 
an angle tracking loop, there are provided means for 
deriving from a target echo an intermediate frequency 
(I. P.) signal, means for estimating an I.F. target 
20 signal frequency periodically, digital filter means 
providing a frequency analysis of signals in the I. P. 
band, the characteristic of the filter means comprising 
a plurality of similar, sequential, overlapping peaking 
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bins, means for prod^se^ng from target signals output in 
adjacent frequency bins in the vicinity of the estimated 
target signal frequency a derived signal characteristic 
having a peak at a predetermined position in relation 
5 to the estimated target signal frequency, and means for 
correcting the estimate of target signal frequency and 
shifting the derived signal characteristic accordingly, 
the target signal resulting from said derived signal 
characteristic being employed in a said tracking loop. 
10 Where there is included a doppler tracking 

loop having a speedgate filter in the I.F. signal path, 
the derived signal characteristic may have sections 
centred symmetrically above and below the estimated 
target signal frequency, means then being provided for 
15 comparing target signal constitutents of the respective 
sections and thereby determining the frequency error 
between the actual and estimated target signal frequency, 
the frequency error being employed to control the I. p. 
target frequency to tend to maintain it at a predeter- 
mined frequency within the pass band of the speedgate 
filter and to tend to bring the derived characteristic 
into alignment with the actual target signal frequency. 

In a missile-borne radar tracking system as 
aforesaid the width of the frequency bins is preferably 
25 controllable in dependence upon the estimated value of 
the difference between the target signal frequency and 
the predetermined frequency within the speedgate pass 
band, the bin width being increased in response to a 
high estimated value to give good tracking ability and 
30 decreased in response to a low estimated value to give 
good velocity discrimination against targets moving at 
low relative velocities. _ 

In a radar tracking system including an angle 
tracking loop employing sum and difference channels and 
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and means responsive to a^rlatio of the sum and difference 
channels and means responsive to a ratio of the sum and 
difference signals to provide an indication of angular 
error between target line of sight and boresight, each 
5 of the sum and difference channels preferably employs 

means for producing a derived signal characteristic, each 
such derived signal characteristic comprising a peaking 
characteristic similar to the individual characteristics 
of the digital filter means and centred on the estimated 

10 target signal frequency so as to provide a narrow 

frequency pass band which tracks the estimated target 
signal frequency* 

According to another aspect of the invention, in 
a monopulse radar tracking system employing sum and 

15 difference signals for the determination of target 

direction, the system comprising a doppler tracking loop 
maintaining an intermediate frequency target signal 
within a speedgate filter pass band, and digital filter 
means in each of the sum and difference channels providing 

20 at periodic update intervals an analysis over a plurality 
of adjacent frequency bins of target signal components 
within the said pass band,ancLlna method of confirming 
the presence or absence of a target signal within a 
particular frequency bin, a series of comparison processes 

25 are performed, each involving the sum of a predetermined 
number of successive power output values from the 
particular frequency bin accumulated with any previous 
such sums and a comparison of such cumulative sum with 
upper and lower threshold values which become progressively 

30 closer with each comparison process, confirmation of the 
presence of a target signal in any comparison process 
being indicated when the upper^hreshold is exceeded by 
the cumulative sum, confirmation 6f the absence of a 

35 target signal being indicated whj$n the lower threshold 
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exceeds the cumulativ^sum, and*' a further comparison 
process being initiated wherf*the cumulative sum lies 
between the upper and Jower thresholds. 

In such a method, a running average may be 
established incorporating a fixed number of bin power 
output values, the earliest incorporated bin power 
output value being discarded as a current value is 
incorporated so as to produce a running average bin 
power, this running average being compared with a 
predetermined threshold to provide a target-signal- 
present or absent indication the cumulative sum 
indication and the running average indication 
contributing to a net condlusion in which the cumulative 
sum indication takes precedence if the running average 
indication points to a target signal absent, and in 
which the possible results of the cumulative sum 
indication: target signal absent, indeterminate, 
present, are treated as: indeterminate, present, present 
respectively if the running average indication points 
to a target signal present. 

According to a further aspect of the invention, 
a radar tracking system employing sum and difference 
signals for the determination of target direction, 
includes digital filter means in each of the sum and 
difference channels providing at periodic update 
intervals an analysis over a plurality of adjacent 
frequency bins of potential target signal components, 
means for identifying a target frequency bin, means for 
applying sum and difference signals in respect of the 
identified target frequency bin as input signals to 
product means for producing a complex product of one 
of the input signals and the complex conjugate of the 
other, a signal-to-noise ra-tffo indication being derived 
from the imaginary component o£$tfte complex product and 
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a power level indication of>tfhe sum channel signal 
within the target £reqq$sncy bin, the signal-to-noise 
ratio indication thus having a high value in the presence 
of incoherent reflections from multiple targets and a low 
5 value in the presence of coherent reflections from a 
single target, the system further including a basic 
signal-to-noise ratio indication derived from power level 
within the target frequency bin and average power level 
over the plurality of frequency bins, this basic 

10 indication thus having a high value in the presence of 
single or multiple targets within the target frequency 
bin and a low value in the presence of wideband noise, 
and means responsive to both of the signal-to-noise 
ratio indications to provide an indication of single 

15 or multiple targets. 

. The relevant f eatures of one embodiment of a 
pulse doppler radar seeker will now be described, by 
way of example, with reference to the accompanying 
drawings, of which ;- 

20 

Figures 1, 2 and 3 together show, in block 
diagram f orm, the basic elements of the system, 

Figure 2 showing a doppler tracking loop, and 

Figure 3 showing an angle tracking loop 
25 together with single/multiple targer discrimination 
circuitry; 

Figure 4 is a table illustrating a combination 
of target confirmation tests; 

and Figure 5 is a block diagram of an angle 
30 tracking and servo-control system. 
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In that part of the tracking system shown in 
Figure 1, an antenna $^ransmits radar pulses of 
frequency controlled by a local oscillator 4- and doppler 
tracking control path 6, the local oscillator and 
5 control signals being combined by a mixer 8 and applied 
to the antenna by way of an amplifier 9 and circulator 
11, in known manner. The antenna includes a so-called 
comparator circuit which takes the target echo (and any 
other received signals) as received by a square array 

10 of four elements of the antenna and produces three 

output channels, the azimuth and elevation difference 
channels 13 and 15 and the sum channel 17 , again in 
known manner. 

The antenna is steerable with respect to the 

15 missile by means of servo-controlled motors 19 to attempt 
to maintain the antenna boresight directed at a target* 
The actual direction of the boresight relative to the 
missile is determined by gimbal •pick-off 1 transducers 
21. 

20 TI}e sum and difference signals are applied to 

respective mixers 25, 25 and 27 which each have an 
input from the local oscillator 4. The transmitted 
frequency is the sum of the local oscillator frequency, 
a basic intermediate (I. P.) frequency derived from a 

25 voltage controlled oscillator (VCO) 29, and the estimated 
doppler frequency f d corresponding to the voltage applied 
to the VCO 29 and being the output of the doppler- 
tracking loop. The output of the mixers 23, 25 and 27 
is therefore the I.F. plus the error in the estimation 

30 of the doppler frequency (true doppler frequency = f d ). 

The mixers 23, 25 and 27 are followed by head 
amplifiers 33, 35 and 37 • * T r 

Electronic angle trackiri^ (E.A.T. ) circuitry 

; CO-NFIBENTIAL 
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then follows in which a^espective proportion of the sum 

signal is added into each of the difference channels, 

the proportion being k for the azimuth difference and 

a 

k for the elevation difference. This proportion is 

5 calculated to produce a net zero difference channel 

signal and thus simulate alignment of the bores ight and 

target, the actual of f-boresight angle information (apart 

from error in the derivation) being contained in the 

control signals £ Q and £ (estimated boresight/target 

sl e 

10 sightline error) applied to the E.A.T V elements 39 and 
41 on lines 47 and 49. 

The resulting sum and difference signals are 
applied to range gates 53 9 55 and 57 which are open for 
a predetermined short time at a controlled delay 

15 period after each radar pulse transmission. The timing 
of this delay period is controlled by a range tracking 
loop circuit 59. The sum channel is branched to feed an 
auxiliary range gate 51 which differs from the others in 
that the gated period is in two halves in one of which 

20 the gated signal is inverted. The output of the sum 
range gate 57 and the part inverted output of the 
auxiliary range gate 51 are applied to a phase-sensitive- 
detector which thus gives a zero output if the received 
signal lies equally in the first and second halves of the 

25 auxiliary gate. If the signal lies earlier or later, 
the PSD output will have a net negative or positive 
value which is used by the RTL circuit 59 to re-centre 
all of the range gates on the target signal. The target 
is thus effectively range tracked. 

30 The range tracking loop 59 also applies an 

inhibit signal to the transmitter amplifier 9 while the 
range gates are open to ensipje that no breakthrough of 
transmitted pulse through iJhe circulator 11 can interfere 



€■ 



\ 



cq.n.fidMtiau^ 

with target echo signals. 

Following the^ange gates 51, 53, 55 and 57, 
^ the sum and difference signals are applied to speedgates 

61, 63, 65 and 67 respectively. These are in fact band 
5 pass filters covering the possible range of doppler 
frequencies of interest and centred on the above 
mentioned intermediate frequency. The doppler tracking 
loop attempts to maintain the various sum and difference 
target frequencies in the centre of the speedgate pass 
10 bands by control of the VCO 29. 

The sum and difference signals are then scaled 
in accordance with the sum channel signal level by 
age amplifiers 71, 73, 75 and 77 controlled by an age 
detector 69. 

15 The signals until this point have all been 

analogue but are now converted by analogue digital 
converters 81, 83, 85 and 87 for the purpose of digital 
frequency analysis. In each of the sum and difference 
channels is a f ast-f ourier- transform (PFT) filter (91, 

20 93, 95 and 97) • Each such filter comprises an array 

of perhaps 32, 64, 128, 256, filter elements, typically 
say, 128, these filter elements respectively covering 
adjacent narrow frequency bands and together covering 
the speedgate passband. The frequency bands of such 

25 an FFT filter are known, and will be referred to in 

this specification, as frequency bins. The characteristic 
of each bin is of peaked form and overlaps with those of 
the adjacent bins. The bin outputs arise in digital 
form, representing the amplitude and phase of the signal 

30 component within the particular bin. One set of such 

output data requires information from 128 analogue samples. 
The output data rate of the FFT filters, i.e. the 
•update rate* is 78 Hertz if the 'PPT. overall bandwidth 
is 10 kHz and there are 128' bins .\ Thus the update rate 

35 is. equal to the individual bin 'w^dtau 
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The output data rate etf^fhe FET>f liters, i.e. the 
•update rate* may thence about 80 Hertz. 

The width of the frequency bins is controllable, 
as indicated by control circuitry 89, thus covering the 
I.F. range with few, broad bins, or a greater number 
of correspondingly narrower bins. 

Each of the sum, difference, and range channels 
is thus analysed into 128 (say) frequency bands within 
the IP doppler range and a target signal frequency can 
be identified very narrowly. 

Referring now to Figure 2, this shows the 
remainder of the doppler tracking loop, which is -entirely 
digital and is constituted by processes performed by a 
data processor. The input is derived on line 97' from 
the sum channel FFT filter 97 in Figure 1, the single 
input 97 • representing all of the 128 sum bin outputs. 
A target bin detection process 101 makes an initial 
estimate of the target bin identity, as will be explained, 
giving an initial target bin frequency f Q , (as a 
displacement from the speedgate centre frequency). 

A sliding filter process 105 to which the FFT 
filter-outputs are applied, produces, in response to the 
initial estimate of target frequency f Q , a derived 
characteristic comprising two individual bin 
characteristics disposed symmetrically above and below 
f Q . One of these bins is derived from the left and 
centre ones of three adjacent bins and the other from 
the centre and right of the three. Discriminator (105), 
age (107) and divider (109) processes produce an error 
signal 6 d , being the frequency error between the estimated 
target frequency (f Q in this case) and the actual target 
frequency. 

The error signal £ d is processed further, as 
will be explained, but in^.ad^fii^on is directly added to 
the current target frequency estMate £ (initially 



f Q ) in a process 111. ^*The r^stalt is therefore the current 
magnitude of the targe t^Ibppler frequency relative to the 
speedgate I.F. centre frequency. This digital value is 
subjected to a gain control process G3 (113) and provided 
5 as a digital output of the doppler tracking loop. 

Referring again to Figure 1, the DTL output 
on line 113' is converted back to analogue form and 
applied to an integrator 121 which accumulates the doppler 
frequency error and controls the oscillator VCO (29) 

10 accordingly. With a constant frequency error, i.e. 

constant target acceleration, the loop will lock on to 
the target frequency and cause the VCO frequence to ramp 
in tracking it. At constant relative target velocity 
the frequency error will cease, the VCO input will be 

15 zero and the VCO output will remain constant at the IF. 

An accelerometer 117 detects acceleration of . 
the missile and adds in a corresponding factor to the 
integrator input by way of adder circuit 119. 

Reverting to Figure 2, the further processing 

20 of the discriminator output, referred to above, is also 
employed to produce from a final addition 323 a signal 
f which is a further estimate of the target frequency 
to speedgate centre frequency error. This signal f is 
employed in a target confirmation process 125 to be 

25 explained. 

Confirmation of the target bin is employed in 
several ways: the electronic angle tracking process 
referred to in Figure 1 is enabled; the guidance system 
of the missile is enabled; and incidental to the main 

30 arrangement a signal-to-noise indication is obtained 

for employment in multiple target detection, as will be 
explained. ^ % 

Considering Figure 2 in\reater detail now. 
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and in particular the target*"bin detection process 101. 
it will be recalled th#b the PPT circuits in each sum 
and difference channel provide outputs giving the complex 
amplitude in each of the PFT filter 'bins 1 (typically 
5 16, 32, 64, 128 or 256 'bins' according to the bin width 
control setting 89). The filter bin width is equal to 
the total FFT bandwidth divided by the number of 
filter bins and the output data rate is equal approx- 
imately to the bin width, 

10 Target detection is carried out on all of the 

power outputs (i.e. the square of the sum of the real 
and imaginary components) of the sum channel PFT with 
the exception of the marginal bins at each end, i.e. 
the first and last N L where N L = 9* 5, 3, 2, 2 for 

15 256, 128, 64, 32, 16 bin FFTs respectively. The power 
output of each bin is divided by the average of all the 
bins excluding the marginal first and last N L , and then 
this value is compared to a threshold (T D ) which is 
determined so as to allow a certain number of crossings 

20 if the input is pure thermal noise (false alarms). The 
technique used is to compare each bin power, scaled by 
the average power, as above, with the maximum scaled value 
of the previous bins and with the threshold T D . The 
maximum bin at the end of the process which also exceeds 

25 T D is taken as a target alarm for the bin in question. 
If no bin exceeds the threshold then no alarm is found 
and the detection process is applied to the next set of 
FFT data. This detection method can be modified in 
practice to allow multiple alarms. 

30 If the target alarm is found* its bin number 

and centre frequency f Q is^jatput so that the doppler 
tracking loop and confarmationVprocesses can act on the 



correct frequencies. An additional calculation is 
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performed in tlie detection process. The contents of 
the alarm bin are reduced by the above threshold value Tp. 
This is done because the probability distribution of 
the maximum bin which exceeds a given threshold is 
5 approximately a constant (T^) plus a Bayleigh 

distribution (provided T D is > log e N). By subtracting 
the threshold value, the initial alarm can be treated 
in the same way as subsequent power outputs from the bin 
and this simplifies the confirmation system. 
10 The calculations performed in the above process 

for each bin power p ± in order to establish the maximum 
power bin are: 

N-N L 




15 where r^ is the scaled power of a bin numbered i; N is 
the total number of bins; N L is the number of marginal 
bins; and j is the number of a bin within those considered 
for averaging. If x ± > T D and is also > r c max then 
r c is set " r i and 1 ls stored as the alarm bin 

20 number. r c max is the maximum value of r ± •so far*. The 
stored value of i which remains when all bins have been 
considered is then the target alarm bin. This initial 
setting, or estimate, of the target frequency (f in 
Figure 2) is then employed to close the doppler tracking 

25 loop and initiate the target confirmation process. It 
will be appreciated that the target frequency, while 
being within the detected bin^(i), will not in general be 
precisely the bin centre ^frequency f • 

In order- to get a high> probability of target 



f0 



CONFIDENTIAL 



-13- 

acquisition with a low probability of 'false confirm 1 
it is necessary to sum incoherently over a large number 
of FFT update periods. To make this possible, target 
tracking is necessary, since either the target must be 
kept in a fixed bin or the alarm bin must be known in 
some other way. For this reason a two stage acquisition 
process is used. The first stage which has already been 
described detects the whereabouts of a likely target, 
the second, confirmation, ensures that this detected 
signal is indeed a target with a high degree of 
certainty. 

The method of operation of the confirmation 
process, referenced 125 in Figure 2, is as follows 
the 'alarm' bin f is designated by the doppler tracking 
algorithm of Figure 2, except immediately after the first 
'detect when the detection subroutine provides this 
information. As in the detection algorithm, the 
equivalent of Equation (1) above is used to provide a 
scaled target/alarm bin output. (If the alarm bin is 
out of the range Nj+l to N-N^ the scaled bin value is 
set to zero). Then this scaled 'target bin' value on 
one FFT data processing is added to the same quantity on 
the successive update. This summing continues until a 
predetermined number of FFT power output values from 
the target bin have been added together. At this point 
the sum is compared to two thresholds, an upper and lower 
threshold. If the sum of the bins is greater than the 
upper threshold then a 'target present' flag is set, 
if the sum is less than the lower threshold, 'target 
absent' is confirmed, or in other words the suspected 
target is rejected. If the sum is in between then the 
flag remains in a 'don't know' state. In this case a 
further set of 'target 1 bin power values from successive 
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FFT updates are added to the first set and the accumulated 
total compared to two new and closer thresholds, the upper 
one again denoting •confirmed presence 1 and the lower 
one 'confirmed absence 1 . 
5 After a predetermined number of such comparison 

stages (provided as an input to the program) the upper 
and lower thresholds are made to coincide so that a 
definite decision is forced although it may hot be 
conclusive, as will be seen. 

10 This successive comparison test is called the 

•cumulative stun 1 indication. 

At this point a running average of the last N 
sets, each of N PPT updates begins. As the current bin 
power value is incorporated in the running total the 

15 earliest value so incorporated is discarded. The 

criterion for confirmation is that the running average 
should be above a threshold (normally taken to be the 
last comparison process converged- threshold) and 
that in addition the individual comparison processes 

20 should not reject. This latter test is put in to ensure 
reasonably rapid response if the target should suddenly 
disappear for some reason. 

The result is a 'running average 1 indication • 
The •cumulative sum 1 indication and the 'running 

25 average 1 indication are combined as illustrated in Figure 
4 to provide a net conclusion. It may be seen that the 
cumulative sum indication takes precedence if the running 
average indication is negative, while the cumulative sum 
indication is, in effect, endorsed by one level of 

30 indication certainty if the running average indication 
is positive. 

In addition a certain amount of re-initialisation - 
takes place. If the overall confirmation goes off the 
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•confirm present 1 state the running average is stopped. 
Also if the cumulative sum test rejects or confirms, the 
accumulation of sub-set sums ceases and each subsequent 
sub-set is tested separately. 
5 When target confirmation (125) is complete the 

•confirm present 1 output 126 is used to enable the angle 
tracking loop and the missile guidance loop. 

The doppler tracking system employed allows 
doppler tracking to take place both within the digital 

10 processor of Figure 2 and through the analogue VCO 29 
in Figure 1. The key function that allows tracking to 
take place within the data processor itself is the 
sliding filter /discriminator process, 103, 105, 107, 109. 
To form a frequency discriminator the power 

15 output of two neighbouring bins can be subtracted* If 

the target signal sits symmetrically astride the junction 
of the two bins the power in each will be the same and 
the difference will be zero. If the target signal is 
off centre one way or the other the result will be 

20 positive or negative accordingly. If two such bins, 
centred on a specified target frequency estimate, can 
be simulated, therefore, frequency shift from this 
central position can be detected and the target signal 
frequency tracked continuously. 

25 To vary the position of the discriminator 

continuously in this way a technique of (effectively) 
sliding the filters of the FFT filter 97 is adopted, 
which requires the simulation of each of the two 
simulated bins mentioned above. The following applies 

30 to the simulation of each of the two. Given the complex 
output of two existing adjacent FFT bins it is possible 
to construct by using these two quantities alone a new 
FFT bin which will have as its peak any chosen value of 
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frequency between the mid-point of the two original bins. 
In order to derive a simple algorithm a rectangular 
^ window is assumed and phase factors between adjacent 

bins of n /B are ignored. The new FFT bin characteristic 
5 can be written approximately as: 

f( X ) = a si ; tc x + p 8i g 1 ; x ^ 1 " x? (2) 

where x is the fraction of a bin from the centre of the 
first FFT bin where the response is required., a is -the - - 
amplitude contribution of the first bin, and 3 is the 
10 amplitude contribution of the second bin. It can be 
shown that if a, p are related by a parameter 6 as 
follows: 



a 



- 6 2 fl + 7tcotTi6(l-6)3 



7icotu6(6(l-6)(l-26)) - 1 + 26 - 26 2 
P = 1 - a 

15 then the point f (x = 6) is a maximum of the filter of 

Equation (2). Hence a sliding filter can be implemented:, 
by taking the frequency f, evaluating the bins nearest to 
this and then calculating the parameters a, p needed to 
give two filter bins each separated by one half a bin 

20 separation from the estimated target frequency (2). 

The same technique is used on the filters regardless of 
the type of window that is being used. The equations 
for the sliding filters thus become :- 

N t = Integral part of ((f + 0.5B)t + 1.5) 

oca S 

C 1 » oc(x) (i.e. the above function a 

calculated for predetermined 
set. of x values) 
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to give optimum performance. The age is thus 
agc Q = agc n-1 k + (1-k) /°agc/ 



where S is another scale correcting factor for the 
c 

different gains found at different positions of the 

5 sliding filter and different windows. This is given by 

— —2 
= 1 + x b n + x b 0 
c -L d 

where b^ and bg are constants depending on the type of 

window. k equals exp (-t e /t ) where t e „ is the age 

s ago age 

time constant. The resultant discriminator output 
10 (from divider 109) is 

£ d * d S c V^n *• > Hz 

The digital tracking loop is closed entirely 

within the data processor (see Figure 2) and provides 

the estimate of target frequency f that is needed for 
15 the confirmation process. The particular implementation 

used is a type two loop. A block diagram of the system 

is shown in Figure 2. 

It can be seen that there are two outputs of the 

loop, one for the confirmation process (f) and one 
20 feedback (f) to the sliding filter/discriminator. The 

transfer functions of these two are different and are 

given by 

£ zGjU+Gg) - Gj^ 

7 z 2 + z(G 1 (l+G 2 )-2) + (l-G^ 

~ z 2 G ± G 1 (G 2 -l)z 

T = ~2 



z c + z(G 1 (l+G 2 )-2) + (1-G X ) 

25 Where f is the target frequency relative to the centre 
of the IP sum channel. These transfer functions are 
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C a - C l ^ta" 1 ) " ^"V b < N ta> 
C b = (l-C^) b(N ta+ l) - C x b(N ta ) 

where lff^ is the number of the bin 'intersected 1 by the 
required frequency, f is the estimated doppler frequency 
5 according to the digital tracking loop, B is the overall 
FFT bandwidth, t Q is the update time(= No, of bins-rB), 
N is half the number of bins in the PPT and b(N) is the 

P 4-V, 

the N PFT bin (complex). Thus C and C. are the 

& D 

sliding filter outputs of simulated adjacent bins centred 
10 on the frequency f , and are again complex. 

It is then necessary to construct the 
discriminator 105 of Figure 2, and this can be done by 
taking the frequency error d « j C a\ 2 " | c t> | 2# 
However, the sliding filter technique results in a 
15 scaling at the origin of this quantity which depends on 
filter position and hence a correction for this has to 
be applied. A quadratic correction term is used of the 
form 



20 



S d " a l + a 2 * + a e * 2 



where x = min (x, 1-x) and a 1$ a 2 , a^ take various values 
for the different window functions being used. This also 
allows the discriminator to be automatically scaled for 
the appropriate window function. 

The discriminator output 'd 1 must also be 
25 automatically gain controlled and this is accomplished 
by a simple first order feedback system operating off a 
square law detector 107* The detector operates off a 
new bin of the form 

age a b 

which gives a filter centred on the target frequency • 
In conjunction with the discriminator this is believed 
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implemented by taking the output of the discriminator t d 
and then multiplying by the gain (Gl). 

The algorithms for the implementation of these 
two closed loop transfer functions are as follows. If 
the discriminator output is £ ^ then 

f^n) = f-^n-l) + G 1 £ d 

f 2 (n) = f 2 (n-l) + G^Cn) 
f (n) = f-^n) + f 2 (n) 

f (n) = f-^n) + f 2 (n-l) 

These algorithms are illustrated in Figure 2 as 
follows. The output £ d of. the discriminator processes 
105, 107, 109 is subjected to a gain Gl (127) producing 
an input G^ £ d to a summing process 129 a storage delay 
device 131 feeds back the output of the summing process 
129 to its input at the following update. If f^n) is 
the output of the summing process 129 on the nth update, 
f^n) must therefore be equal to G x £ d + f^n-1), the 
loop therefore constituting a digital integration. 
The output of the loop, f^n) is applied to a further 
storage delay device 133 whose output f^n-1) is applied 
to a summing process 135* A second input to this summing 
process is derived from the output f^Cn) of the above 
loop as subjected to a gain G2 (137) to give Ggf 1 (n) and 
a further integrating loop 139, 14-1 to give f 2 (n) equal 
to f 2 (n-l) + G 2 f 1 (n). A further storage/delay device 143 
then gives f 2 (n-l) the other input to summing process 135. 

The output of summing process 135 is thus 
f(n-l) = f^n-1) + f 2 (n-l) i.e. the estimated target 
frequency based on the preceding update of the FFT data. 
This output is applied to the sliding filter process 103 
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for calculation of the next error and is also 

d 

applied to summing process 111 together with the current 
error £ d as previously described, to provide the doppler 
tracking loop output. 

The signal f 2 (n-l) is applied to a further 
summing process 123 hut in this case with the current 
input to process 135, i.e. f-^n). The summed output of 
process 123 is thus f^n) + f 2 (n-l), denoted f and 
applied as the estimated target frequency to confirmation 
process 125, as previously described. 

Referring now to Figure 3 the FFT outputs from 
the azimuth difference filters 95 are input on lines 
represented by 95* and the sum FFT outputs similarly 
on lines 97 These sum and difference outputs, each 
from, typically, 128 bins, are applied to respective 
sliding filter processes 145 and 147 similar to those 
described for the doppler loop of Figure 2. In each 
case the estimated target frequency f derived from the 
doppler loop is used to select the bin in which it falls, 
a new or 'supposed 1 bin being derived from this bin and 
an adjacent bin, the simulated bin having its centre or 
peak frequency aligned with the estimated target frequency. 
Changes in the target frequency as determined by the 
doppler tracking loop cause the value of f to change, 
the simulated bin remaining locked to it and thus sliding 
up and down the I.F. range with the doppler error 
frequency and providing a high degree of target velocity 
discrimination . 

The sliding filter outputs are in digital form 
and are therefore applied to the digital equivalent of 
a phase-sensitive-detector. Thus the difference signal 
in complex form has its complex conjugate formed by a 
process 149. The sum and complex conjugate of the 
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diff erence are then multiplied together by a process 151 
which produces real and imaginary outputs. The real part 
is divided by an age detector (157) output to give the 
amplitude ratio for each of the two channels i.e. 

(D e -k e S)S (D a -* a S)S 
j and 

where S , the mean square of the sum signal, is the output 
of the age detector 157. Only the azimuth difference 
channel is shown but the elevation difference channel is 
treated similarly. This is fed to a Kalman filter in 

10 Figure 5 to give an output £ • £ is added to an 

aberration correction factor and fed to the EAT gain 
element to close the loop. The Kalman filter in Figure 
5 produces an estimate £ a of the boresight error, that is, 
the estimated value of the target angle off boresight. 

15 This error estimate is combined in a summing process 
156 with a boresight direction indication relative to 
the missile. The latter is provided basically by gimbal 
pick-off transducers 21 (shown in Figure 1 also) 
corrected for radome aberration (150) which itself is 

20 dependent upon; boresight direction relative to the 

radome; radome temperature (163); and transmitted 

frequency (161). The aberration corrector thus produces 

azimuth and elevation angles of the effective boresight 

rather than the physical boresight. 

25 The resulting signal k output from summing 

a 0 

process 156 is employed in Figure 1 to control the EAT 
element 41 and close the EAT loop. 
A When the EAT loop is in equilibrium the value 

6 gives a true measure of the boresight error. 
30 Considering now the imaginary output of the 
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product 151, this is divided by the age signal in divider 
155 to give a form of signal-to-noise indication. This 
output of the divider 153 gives an indication of the 
power of the signal components that are either incoherent 
between the sum and difference channels or in phase 
quadrature. The 'incoherent' power component is large 
when thermal noise or jamming in either the sum or the 
difference channels is large, and tends to zero at high 
signal to noise ratio. The incoherent power output is 
also large on extended targets or on targets flying in 
formation that are unresolved by the seeker. 

The signal-to-noise (S/N) signal output from 
153 is applied to a dividing process 159 the other 
input of which is derived from the S/N output of the 
15 target confirmation process 125 of Figures 2 and 3.. 
The two S/N indications differ in that the basic 
indication, from the confirmation circuit is derived 
as a ratio of the power in the target bin to the average 
power over the FPT band. It does not therefore 
20 distinguish between a coherent single target in the 
target bin and incoherent multiple targets within the 
target bin. TheQPSD S/N indication is derived solely 
from the target bin and distinguishes only between 
coherent (single) targets on one hand and multiple 
25 targets and noise on the other, giving a low value for 
the coherent target and a high value for the incoherent 
targets. If, therefore, the qpSD S/N indication is 
divided by the basic S/N indication from the confirmation 
process the result will be a low value for a single 
30 coherent target in the target bin and a high value for 
multiple targets in the target bin. The result of the 
division process 159 thus gives an indication of single 
or multiple targets. 
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The system as described is adaptable for the 
various missile/target conditions that may arise in an 
engagement, and in particular in the terminal stages 
where conditions may change rapidly. 

The parameters in the doppler loop are made 
adaptive to allow for changing target conditions so 
that optimum performance is maintained at all times. 
Referring to Figure 2: 

(i) Immediately after target detection, gains Gl, 
G2 and G3 are set to high values so that the loop 
settles quickly. The gains are then reduced until, 
after a predetermined length of time, the gains take 
constant lower values to reduce the noise in the loop. 
The gains are increased as the missile flight progresses 
according to the signal strength measured in the age 
amplifiers (connection from age not shown), or by one 
or other of the signal to noise assessment indications 
described above, and according to the FFT bin width 
selected (see below). 

(ii) During the missile boost phase, the doppler 
tracking loop experiences substantial transients and the 
FFT bin width is selected to be a large value. Gains 
Gl, G2 and G3 are also set to high values to give a 
rapid response. 

(iii) If a target at long range is to be acquired, 
the FFT bin width is set to be very narrow, and Gl, G2 
and G3 are set to low values after loop pull in. This 
optimises the acquisition threshold. 

(iv) If the target is at short range the FFT bin- 
width is set to be wide and Gl, G2 and G3 are set to 
high values. This optimises acquisition speed and 
minimises tracking errors, A wide FFT binwidth gives 
a faster update rate and hence gives the optimum 
performance very close to impact when rapid response is 
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important. Target range is determined from the range 
tracking loop (Figure 1). 

(v) Once the loop has pulled-in, the signal f is 
proportional to the target frequency rate, i.e. the 
5 target acceleration. If f is large the FFT binwidth is 
set to a high value to allow a rapid tracking ability. 
The loss in doppler discrimination against other targets 
at similar speeds is not important since a manoeuvring 
target is not able to keep in close proximity to other 

10 targets. When f is small, a small FFT binwidth is used 
to give maximum doppler discrimination against targets 
flying in formation. 

(vi) When multiple targets in close proximity are 
detected (e.g. by the output from divider 159 in Fig. 3) 

15 or the target acceleration is small , the FFT binwidth 
is set to a low value to give maximum doppler 
discrimination. 

The overall angle tracking system can be 
conveniently divided into three parts. Firstly the EAT 

20 receiver, secondly the angle tracking filter and thirdly 
the servo control and stabilisation system. 

A block diagram of the complete angle tracking 
and servo control system is shown in Figure 5. The 
angle tracking filter estimate c of boresight error 

25 angle in each of the azimuth and elevation channels is 
used to add a controlled portion of the sum channel 
which corresponds to these angles according to the 
stored D/S slope at the origin, into the difference 
channels, as described above with reference to Figures 

30 1 and 3* These difference signals then pass through 
similar IF chains as for the sum channel, are sampled 
and analogue-to-digital converted (81-87) • T hey are 
then fourier transformed (91-97) to give the spectrum 
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of signals on the azimuth and elevation channels. The 
sum and difference channels are divided by age signals 
digitally and elevation and azimuth phase-sensitive- 
detector outputs are formed (155) by taking the real 
part of D*xS; these outputs are proportional to the 
angular error between the true target boresight error 
and the best estimate of the boresight error. These 
calculations are all performed at the same update rate 
as the doppler and range tracking loops (i.e. at the 
update rate of the FFT output). 

After the PSD's (151) it is necessary to 
perform an axis transformation in order to transform 
from receiver axes to line of sight axes. This is 
implemented digitally. These signals are then passed 
through an angle tracking (Kalman) filter (171) which 
produces best estimates of sight line rate using a 
variable gain algorithm. The sight line rate signals 
are used in various ways; firstly they provide an 
output to the autopilot as the guidance command, for 
which purpose they are transformed into missile axes; 
they are also used as the input to the EAT integrator 
which feeds estimate of angle back to the EAT element 
and finally it is used as input to the line of sight 
observer (173) to drive the antenna mechanism. This 
latter is added to an inert ial reference unit (IBU) 
gyro output (175) in order to provide the total movement 
of the target as measured in line of sight axes. 

These angular rates are still calculated at 
the (low) update rate of the loop. These rates are 
t£en integrated using a high update digital integrator 
177 and the output is transformed from line of sight 
axes into motor shaft angles using digital axes 
transformations 179* The angle demands are then used 
to control a position control loop 181 for the antenna. 
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The angle demand 183 is subtracted from the potentiometer 
pick-offs 21 and the error is fed into a controller 185 
which feeds demands to the antenna mechanism, servo and 
motors 187* The antenna/reflector plate is thus position 
slaved to the integrated IRU gyro and the target sight 
line error. A high update rate controller is used 
The effect of this is to remove body motions on the 
receiver output. An additional line of sight feed back 
189 is used to the EAT element (39, 41), This is added 
to the other feedback path and ensures that the best 
estimate of sight line error can be used to subtract 
from the actual sight line error for EAT, The combined 
EAT feedback path has to be angle transformed (191) 
back into receiver axes and a radome aberration 
correction (159) is added in, as a function of the 
gimbal angles , as described in relation to Figure 3. 

The Kalman filter parameters are varied as a 
function of conditions during missile flight and in 
general are arranged to give high bandwidth when the 
estimated time to impact is small* Time to impact is 
estimated by dividing the missile-target range (from 
the range tracking loop) by the closing velocity (from 
the doppler tracking loop). * The filter bandwidth is 
also made a function of signal to noise ratio, the band- 
width being low when the signal to noise ratio is low. 

It may be seen that a versatile tracking system 
has been described. The sliding filters permit 
discrimination between closely spaced targets the doppler 
resolution being twice as good as with fixed FFT filters. 

An improved signal/noise ratio can also be 
obtained by sliding the simulated bin peak on to the 
target. 

Discrimination of multiple targets, e.g. by 
tracking one target and monitoring another, can be 
achieved merely by duplicating the sliding filters. 
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CLAIMS 

!• A monopulse radar tracking system, including 

either or both of a doppler tracking loop and an angle 
tracking loop, means for deriving from a target echo 
an intermediate frequency (I. P. ) signal, means for 
estimating an I. P. target signal frequency periodically, 
digital filter means providing a frequency analysis of 
signals in the I.F. hand, the characteristic of said 
filter means comprising a plurality of similar, sequential, 
overlapping > peaking characteristics, defining respective 
adjoining frequency bins, and means for producing from 
target signals output in adjacent frequency bins in the 
vicinity of the estimated target signal frequency a 
derived signal characteristic having a peak at a 
predetermined position in relation to said estimated 
target signal frequency, means for correcting the 
estimate of target signal frequency and shifting said 
derived signal characteristic accordingly, the target 
signal resulting from said derived signal characteristic 
being employed in a said tracking loop* 
2* A radar tracking system according to Claim 1, 

including a doppler tracking, loop having a speedgate 
filter in the I.F. signal path, wherein said derived 
signal characteristic has sections centred symmetrically 
above and below the estimated target signal frequency, 
and means are provided for comparing target signal 
constitutents of the respective sections and thereby 
determining the frequency error between the actual and 
estimated target signal frequency, said frequency error 
being employed to control the I.F. target frequency to 
tend to maintain it at a predetermined frequency within 
the passband of said speedgate filter and to tend to 
bring said derived characteristic into alignment with 
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the actual target signal frequency, 

3* A radar tracking system according to Claim 2, 

wherein said derived characteristic comprises two, 
overlapping, peaking characteristics similar to the 
individual characteristics of said digital filter means, 
the two bin characteristics thus derived having peaks 
displaced on opposite sides of the estimated target 
signal frequency by half of a frequency bin width. 

4. A radar tracking system according to Claim 5, 
wherein said means for comparing target signal constituents 
of the respective sections comprises means responsive to 
the power outputs within the two derived bin character- 
istics and means for obtaining an algebraic difference 

of the power outputs to indicate the degree and direction 
of the frequency error* 

5. A radar tracking system according to any of 
Claims 2, 3 and 4-, wherein said frequency error is added 
to the estimated target frequency to provide a 
controlling signal tending to maintain the I. P. target 
signal in the centre of the speedgate pass band. 

A radar tracking system according to any of 
Claims 2-5, wherein said digital filter means comprises 
an array of fast fourier transform filters each providing 
one said frequency bin, the output data being up-dated 
periodically and at each such up-date providing a new 
estimate of the target signal frequency derived from the 
preceding estimate and the current error. 
7« A missile-borne radar tracking system according 

to any of Claims 2 - 6 including means responsive to 
high acceleration rates of the missile to increase the 
width of said frequency bins and thus improve the quality 
of tracking at the expense of target discrimination. 
8 * A radar tracking system according to Claim 7 
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wherein said doppler tracking loop incorporates 
controllable gain elements which are set to high gain 
values in response to high acceleration rates. 
9- A missile-borne radar tracking system according 

to any of Claims 2-8, wherein the width of said frequency 
bins is controllable in dependence upon target range, the 
bin width being set to a low value in response to a long 
range target to optimise the acquisition threshold. 

10. A radar tracking system according to Claim 9, 
wherein said doppler tracking loop incorporates 
controllable gain elements which are set to low gain 
values in response to a long range target. 

11. A missile-borne radar tracking system according 
to any of Claims 2-10 wherein the width of said 
frequency bins is controllable in dependence upon the 
estimated value of the difference between the target 
signal frequency and said predetermined frequency within 
the speedgate pass band, the bin width being increased 
in response to a high estimated value to give good 
tracking ability and decreased in response to a low 
estimated value to give good velocity discrimination 
against targets moving at low relative velocities. 

12. A radar tracking system according to Claim 1 
including an angle tracking loop employing sum and 
difference channels and means responsive to a ratio of 
the sum and difference signals to provide an indication 
of angular error between target line of sight and 
boresight, wherein each of said sum and difference 
channels employs means for producing a said derived 
signal characteristic, each such derived signal 
characteristic comprising a peaking characteristic similar 
to the individual characteristics of said digital filter 
means and centred on the estimated target signal 
frequency so as to provide a narrow frequency pass band 
which tracks the estimated target signal frequency. 
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13 • In a monopulse radar tracking system employing 

sum and difference signals for the determination of 
target direction, the system comprising a doppler 
tracking loop maintaining an intermediate frequency 
target signal within a speedgate filter pass band, and 
digital filter means in each of the sum and difference 
channels providing at periodic update intervals an 
analysis over a plurality of adjacent frequency bins of 
target signal components within the said pass band, a 
method of confirming the presence or absence of a target 
signal within a particular frequency bin, in which a 
series of comparison processes are, performed, each 
involving the sum of a predetermined number of 
successive power output values from the particular 
frequency bin accumulated with any previous such sums 
and a comparison of such cumulative sum with upper and 
lower threshold values which become progressively 
closer with each comparison process, confirmation of the 
presence of a target signal in any comparison process 
being indicated when the upper threshold is exceeded by 
said cumulative sum, confirmation of the absence of a 
target signal being indicated when the lower threshold 
exceeds the cumulative sum, and a further comparison 
process being initiated when the cumulative sum lies 
between the upper and lower thresholds. 
14. A method according to that of Claim 13 wherein 

after a said series of a predetermined number of said 
comparison processes the upper and lower thresholds 
are t .'made coincident so as to force a target-signal- 
present, or absent, indication. 

15« A method according to Claim 13, wherein a 

running average is established incorporating a fixed 
number of said bin power output values, the earliest 
incorporated bin power output value being discarded as a 
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current value is incorporated so as to produce a running 
average bin power, this running average being compared 
with a predetermined threshold to provide a target- 
signal-present, or absent, indication, the cumulative 
sum indication and the running average indication 
contributing to a net conclusion in which the cumulative 
siim indication takes precedence if the running average 
indication points to a target signal absent, and in which 
the possible results of the .cumulative sum indication: 
target signal absent, indeterminate, present, are treated 
as: indeterminate, present, present, respectively if the 
running average indication points to a target signal present 
16. A radar tracking system employing sum and 

difference signals for the determination of target 
direction, the system including digital filter means in 
each of the sum and difference channels providing at 
periodic update intervals an analysis over a plurality 
of adjacent frequency bins of potential target signal 
components, means for identifying a target frequency bin, 
means for applying sum and difference signals in respect 
of the identified target frequency bin as input signals 
to product means for producing a complex product of one 
of said input signals and the complex conjugate of the 
other, a signal-to-noise ratio indication being derived 
from the imaginary component of said complex product and 
a power level indication of the sum channel signal within 
the target frequency bin, said signal-to-noise ratio 
indication thus having a high value in the presence of 
incoherent reflections from multiple targets and a low 
value in the presence of coherent reflections from a 
single target, the system further including a basic 
signal-to-noise ratio indication derived from power level 
within the target frequency bin and average power level 
over the plurality of frequency bins, this basic 
indication thus having a high value in the presence of 
single or multiple targets within the target frequency 
bin and a low value in the presence of wideband noise, 
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and means responsive to both of said signal-to-noise 
ratio indications to provide an indication of single 
or multiple targets. 

17* A monopulse radar tracking system incorporating 

f ast-f ourier-transf orm digital filters providing echo 
signal analysis within a plurality of adjacent frequency 
bins, and means for simulating a frequency bin at a 
selected intermediate point in the range of frequency 
bins, substantially as hereinbefore described with 
reference to the accompanying drawings. 

18. A monopulse radar tracking system incorporating 

means for detecting and confirming the presence of a 
target, substantially as hereinbefore described with 
reference to the accompanying drawings. 

19» A monopulse radar tracking sy6iem incorporating 

means for distinguishing between single and multiple 
targets, substantially as hereinbefore described with 
reference to the accompanying drawings • 
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